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Abstract 

We prove existence of weak solutions for a diffuse interface model for the flow of two 
viscous incompressible Newtonian fluids in a bounded domain in two and three space dimen- 
sions. In contrast to previous works, we study a new model recently developed by Abels, 
Garcke, and Griin for fluids with different densities, which leads to a solenoidal velocity field. 
The model is given by a non-homogeneous Navier-Stokes system with a modified convective 
term coupled to a Cahn-Hilliard system. The density of the mixture depends on an order 
parameter. 
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1 Introduction 

A fundamental problem in fluid dynamics involves changes in topology of interfaces between 
immiscible or partially miscible fluids. Topological transitions such as pinch off and reconnection 
of fluid interfaces are important features of many systems and strongly affect the flow. Classical 
models based on sharp interface approaches typically fail to describe these phenomena. In 
recent years, diffuse interface models turned out to be a promising approach to describe such 
phenomena. In this approach, an order parameter is introduced which allows for a mixing in an 
interfacial zone. Therefore the sharp interface is replaced by a thin interfacial layer in which the 
order parameter, which can be a concentration, rapidly changes its value. Often these diffuse 
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interface models (also called phase field models) allow for local entropy or free energy inequalities 
and in such situations they can be called thermodynamically consistent. It can be justified in 
some cases that sharp interface models are recovered in the limit when the interfacial thickness 
goes to zero. We refer to Lowengrub and Truskinovsky |LT98j . Abels, Garcke and Grtin [AGGllj 
for a result using formally matched asymptotic expansions and Abels and Roger |AR09) for a 
first analytic result. The diffuse interface approach is hence an attractive approach to model 
and to numerically simulate fluid interfaces. 

In the literature, diffuse interface models are well established for two-phase flows of liq- 
uids with identical ("matched") densities, see Hohenberg and Halperin or Gurtin et al. [HH77t 
IGPV96] . In contrast, when the densities are different, several approaches have been discussed in 
the literature. Lowengrub and Truskinovsky | LT98j derived quasi-incompressible models, where 
the corresponding velocity field is not divergence free. On the other hand. Ding et al. |DSS07] 
proposed a model with solenoidal fluid velocities which is not known to be thermodynamically 
consistent. Only recently Abels, Garcke and Griin, see [AGGllj derived a thermodynamically 
consistent diffuse interface model for two phase flow with different densities. It is the goal of 
this paper to show existence of weak solutions for this new model. 

More precisely, we consider the following system of Navier-Stokes/Cahn-Hilliard type: 

dt{p\r) + div(v (g) (pv J)) - dw{2r]{ip)Dv) +Vp= - dw{a{ip)Vip (g) V(/?) in Q, 

div V = in Q, 

dtif + V • Vip = div (m((^)V/u) in Q, 

= ^'{if) + a'(v9)^^^ - div {a{ip)Vip) in Q, 

where J = — f^~P'^ m(y?)V/x, Q = $7 x (0, oo), and C M"', d = 2,3, is a sufficiently smooth 
bounded domain. We close the system with the boundary and initial conditions 

v|af7 = on dfl X (0,oo), 

dn^lan = dnf^ldn = on x (0, oo), 

(v,(^) |4=o = (vo,v^o) in il, 

where dn = n ■ V and n denotes the exterior normal at dQ. Here v and p = p{(p) are the 
(mean) velocity and the density of the mixture of the two fluids, p is the pressure, if is an 
order parameter related to the concentrations of the two fluids, and fi is the chemical potential 
associated to ip. Moreover, Dv = i(Vv-|- Vv-^), r]{ip) > is a viscosity coefficient, and m{ip) > 
is a (non-degenerate) mobility coefficient. Furthermore, ^{(p) is the homogeneous free energy 
density for the mixture and the (total) free energy of the system is given by 

^ ^a(^)^ + vl;(^)j dx, 

for some positive coefficient a{ip). 

An important aspect of our contribution is that we consider a class of singular free energies, 
which includes the homogeneous free energy of the so-called regular solution models used by 
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Cahn and Hilliard [CH58]: 

M/((^) = ^((l + (^)ln(l + <^) + (l-V9)ln(l-<^))-|(^2^ (1.1) 

where Q < 9 < Oc- Mathematically, these singular free energies ensure that the order parameter 
stays in the physically reasonable interval, which is [—1,1] if is the difference of volume 
fractions of both fluids. But this leads to singular terms in the equation for the chemical 
potential. In order to deal with these terms we apply techniques, which were developed in Abels 
and Wilke |AW07j and applied in Abels |Abe09b) to the model for matched densities. 

In the case of matched densities, i.e., = /52, J = 0, the model reduces to the well-known 
model of matched densities discussed in Hohenberg and Halperin [HH77] . In this case existence of 
weak solutions and well-posedness were obtained by Starovoitov [Sta97] . Boyer |Boy99| , Liu and 
Shen jLS03j . and Abels |Abe09bj . Moreover, Boyer |Boy01| considered a diffuse interface model 
for fluids with non-matched densities, which is different from our model. He proved existence 
of strong solutions, locally in time, and existence of global weak solutions if the densities of the 
fluids are sufficiently close. First analytic results for the model by Lowengrub and Truskinovsky 
were obtained in |Abe09a[ lAbellj . For more information on diffuse interface models for two 
phase flows of incompressible fluids we refer to Abels, Garcke and Griin |AGG11| . 

The structure of the article is as follows: In Section [2] we summarize some notation and 
preliminary results. Then, in Section[3l we reformulate our system suitably, define weak solutions 
and state our main result on existence of weak solutions. In Section U] we approximate our 
system with the aid of an implicit time discretization and prove existence of solutions for the 
latter system with the help of the Leray-Schauder principle. Afterwards we pass to the limit in 
Section O and prove our main result on existence for the Cahn-Hilliard/Navier-Stokes system. 
Here in particular the compactness of the velocity and the attainment of the initial data for 
the velocity are non-standard. Finally, in the appendix we discuss the corresponding results on 
existence of weak solutions for a similar model, which was derived in Abels, Garcke and Griin 
[AGGlO] before. 

2 Preliminaries 

We denote a®h = {aibj)fj^^ for a,beR'^ and A^y^ = ^{^ + ^^) for a matrix A G M'^^"'. If X 
is a Banach space and X' is its dual, then 

{f,9)^{f,9)x',x = f{9), feX',geX, 

denotes the duality product. We write X y if X is compactly embedded into Y. Moreover, 
if is a Hilbert space, (• , •)// denotes its inner product. Moreover, we use the abbreviation 

(. , .)m = (• , ■)l'2{A.I)- 

Function spaces: If M C R'^ is measurable, L'^{M), 1 < g < oo denotes the usual Lebesgue- 
space and \\.\\q its norm. Moreover, L''{M;X) denotes the set of all strongly measurable q- 
integrable functions/essentially bounded functions, where X is a Banach space. If M = (a, 6), 
we write for simplicity L''{a,b;X) and L''{a,b). Furthermore, / G Lj^q^([0, oo); X) if and only if 
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/ G L''{0,T; X) for every T > 0. Moreover, L'^^^^{[0,oo); X) denotes the uniformly local variant 
of L'^{0,oo; X) consisting of all strongly measurable /: [0,oo) — t- X such that 

II/IIl«j^^([0,oo);X) = sup \\f\\Li{t,t+l;X) < OO. 

If T < oo, we set Ll^([0,T);X) := L'?(0,r;X). 

Recall that, if X is a Banach space with the Radon-Nikodym property, then 

L'?(M; Xy = L'i\M; X') for every 1 < g < oo 

by means of the duality product {f,g) = Jj^j{f{x), g{x))x\xdx for / G L'^' {M; X'), g G 
L^{M]X). If X is reflexive or X' is separable, then X has the Radon-Nikodym property, 
cf. Diestel and Uhl |DU77] . 

Moreover, we recall the Lemma of Aubin-Lions: If Xq Xi X2 are Banach spaces, 

l<p<oo, 1 < q < CO, and / C M is a bounded interval, then 

G LP(/;Xo) : ^ G L^I-X^)^ L^il-X,). (2.1) 

See J.-L. Lions |Lio69) for the case q > 1 and Simon [Sim87j or Roubfcek [Rou90] for g = 1. 

Let C M'^ be a domain. Then W^{i}), mG No, l<g<oo, denotes the usual L'?-Sobolev 
space, W^oin) the closure of C^in) in WJ^^Q), W-'^iQ) = (Wg?o(f^))', and W-Jp{n) = 
{W^{Q,)y . The L^-Bessel potential spaces are denoted by if*(0), s G M, which are defined by 
restriction of distributions in i?^(M'^) to Q, cf. Triebel |Tri781 Section 4.2.1]. We note that, if 
f2 C M'' is a bounded domain with C^'^-boundary, then there is an extension operator Eq which 
is a bounded hnear operator Eq : W^{Q) — )• W"™(M'^), 1 < p < 00 for ah m G N and -Eq/Ih = / 
for all / G W^{Q,), cf. Stein [St701 Chapter VI, Section 3.2]. This extension operator extends 
to Eq : i?*(Q) H-'iW^), which shows that i7^(J^) is a retract of H^'iW^). Therefore all results 
on interpolation spaces of H'{R'^) carry over to H'iQ). 

Given / G L^(0), we denote by /q = -py f{x) dx its mean value. Moreover, for m G M we 

set 

LJ^)(17) :={/GL''(17) :/n=m}, 1 < g < 00. 
Then for / G Lp'iyt) we observe that 

P^f:=f-fn = f-^^jj{x)dx 

is the orthogonal projection onto L^Qj(r2). Furthermore, we define 

i/Jo) = Hl^^i^) = H\n)nLl,^in), (c,d)^i^^(f,) := (Vc, Vd)^2(f,). 
Then Hj-Q^{Q,) is a Hilbert space due to Poincare's inequality. 

Spaces of solenoidal vector-fields: For a bounded domain C M'^ we denote by C^^{Q) in 
the following the space of all divergence free vector fields in Cq°{Q,)'^ and i^(f^) is its closure in 
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the L^-norm. The corresponding Helmholtz projection is denoted by P^, cf. e.g. Sohr jSoh01| . 
We note that P^jf = f — Vp, where p G W2{^) H L'^q^^O.) is the solution of the weak Neumann 
problem 

{Vp,Vip)n = {f,Vip) for all e (2.2) 

Spaces of continuous vector-fields: In the following let / = [0, T] with < T < oo 
or let / = [0,cxd) if T = oo and let X be a Banach space. Then BC{I;X) is the Ba- 
nach space of all bounded and continuous f: I ^ X equipped with the supremum norm and 
BUC{I; X) is the subspace of all bounded and uniformly continuous functions. Moreover, 
we define BCw{I]X) as the topological vector space of all bounded and weakly continuous 
functions /: / — ?• X. By Cg°(0, T;X) we denote the vector space of all smooth functions 
f:{0,T) X with supp/ CC (0,T). Finally, / e W^{0,T;X), 1 < p < oo, if and only 

if /, ^ G LP{0,T; X), where ^ denotes the vector- valued distributional derivative of /. Fur- 
thermore, Wp ^jq^([0, oo); X) is defined by replacing LP{0,T; X) by L^j^^([0, oo); X) and we set 
H^{0,T;X) = WUO,T;X) and //^^^([O, oo); X) := Wl^,^^{[0,^y, X). Finally, we note: 

Lemma 2.1. Let X,Y be two Banach spaces such that Y ^ X and X' ^ Y' densely. Then 

L°°{r,Y)nBuc{r,x) ^ BC^{r,Y). 

For a proof, see e.g. Abels |Abe09a j . 

Embedding results for interpolation couples: Let {Xq,Xi) be a compatible couple of 
Banach spaces, i.e., there is a Hausdorff topological vector space Z such that Xq,Xi ^ Z, 
cf. Bergh and Lofstrom |BL76j . and let (. , .)[5i] and (. , ■)g^r^ ^^[0, l],l<r<oo, denote the 
complex and real interpolation functor, respectively. If in addition Xi Xq densely, then for 
all 1 < p < oo 

W^{0,T-Xo)nLP{0,T;Xi) ^ BUC{r,{Xo,Xi\_,) (2.3) 

continuously, cf. Amann [Ama951 Chapter III, Theorem 4.10.2]. From this one immediately 
gets 

I^pi,l„,([0,r);Xo) nLl^([o,r);Xi) ^ BUC{P, {Xo,X,\^,J . (2.4) 

A result related to energy inequalities: The following lemma will be useful for passing to 
the limit in energy inequalities. 

Lemma 2.2. Let E: [0,T) — [0,oo), < T < oo, he a lower semi-continuous function and let 
D : (0, T) —7- [0, oo) he an integrahle function. Then 

E{0)ip{0)+[ E{t)ip'{t)dt> f D{t)ip{t)dt (2.5) 
Jo Jo 

holds for all (p G Wl{Q,T) with ip(T) = and if > if and only if 

E{t) + / D{T)dT < E{s) (2.6) 



holds for all s < t < T and almost all < s < T including s = 0. 
For a proof, we refer to Abels |Abe09a| . 
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3 Existence of weak solutions 



In this section we prove an existence result for weak solutions of the following Navier-Stokes/ 
Cahn-Hilliard system for a situation with different densities. The complete system is given by 

dt{pv) + div(/9v (g) v) - div{2r]{ip)Dv) +Vp + div(v (g) J) 





= — div{a{ip)Vip (g) Vcp) 


in 


Q, 


(3.1) 


div V 


= 


in 


Q, 


(3.2) 


dt(f + V • Vip 


= div {m{ip)'Vfi) 


in 


Q, 


(3.3) 




= ^'{^) + a'(<^)^ - div {a{^)Vip) 


in 


Q, 


(3.4) 




= 


on 


s, 


(3.5) 


On flan = dnnlan 


= 


on 


s, 


(3.6) 


(v,v?) |t=o 


= (vcV'o) 


in 




(3.7) 



where J = — ^^^^^ 771(93) V/i. In the above formulation and in the following, we use the abbrevi- 
ations for space-time cylinders Q(s,t) = ^ >^ Qt = Q(o,t) Q = Q(o,oo) analogously 
for the boundary S(^s,t) = di^ ^ is,t), St = 5'(o.t) and S = S(q co)- Equation ()3.5p is the no-slip 
boundary condition for viscous fluids, n is the exterior unit normal on dO,, dnn\dn = means 
that there is no mass flux of the components through the boundary, and dn<p\dn = describes 
a contact angle of 7r/2 of the diffused interface and the boundary of the domain. 
We reformulate the first line suitably. To this end, we first calculate 



■ div{a{ip)Vip (g Vip) = — div {a{(p)V(p)V if — a{ip)V 



Multiplying (fO|) with Vip using '^'{ip)Vip = V (^(93)) and a'{ip)Vip = V (0(9?)) then leads to 
the following identity 

- diY{a{^)V^)V^ = fiV^ - V{^{ip)) - V (a((^)) , 

which gives 



- div{a{(p)V(f (g) = fiVip - V('I'(99)) - V \^a{(f) 

With a new pressure g = p + ^(v?) + a{ip) ^'^p we replace ()3.ip by 

dtipv) + div(/9v (g v) - div(2?7((/?)i:»v) + V^f + div(v J) = fiVip . (3.1') 
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3.1 Assumptions and definition of weak solutions 

In the following we summarize the assumptions needed to formulate the notion of a weak solution 
of (I3.1I1-(I3.7I) and an existence result. 



Assumption 3.1. We assume that n CW', d = 2, 3, is a bounded domain with smooth boundary 
and additionally we impose the following conditions. 

(i) The constitutive relation between density and phase field is given by p{(p) = ^{pi + P2) + 
^{P2 — Pi)^ as derived in Abels, Garcke and Griin \ AGO 11^ . where Pi > are the specific 
constant mass densities of the unmixed fluids and if is the difference of the volume fractions 
of the fluids. 

(ii) We assume a,m £ C^(M), rj G C'^(M) and < mo < a(s), m(s), 17(5) < K for given 
constants mo,K > 0. 

{iii) For the homogeneous free energy density ^' we assume that ^ G C([— 1, 1]) H C^((— 1, 1)) 
such that 

lim ^'(s) = -00 , lim ^'(s) = 00 , ^"(s) > -k for some k G M . (3.8) 

s— 1 s—^l 

[iv] Additionally we impose the condition that lims_).-|-i = +00. 

Remark 3.2. (i) An example for \I' is given by 

^{s) := I ((1 + s) ln(l + s) + {l-s) ln(l - s)) - s€ [-1, 1], (3.9) 

where a,b > 0. We note that ^ is convex if and only if a > b. 

(a) The Assumption \3. 1\ (iv) is needed to reformulate the model to apply results from the theory 
of subdifferentials. 

{Hi) As the solution will lie in the interval [—1,1], we only need the functions a,m,r] on 
this interval. We then extend the functions a, m, rj to the whole of M such that (ii) in 
Assumption \3.1\ is fulfilled. 

Now we can define a weak solution of problem ()3.ip - ()3.7p . 

Definition 3.3. Let T £ (0, 00] and set either I = [0, 00) if T = 00 or I = [0,T] if T < 00, 
vq G L'^{^) and ifo G H^{Q,) with |(/?o| < 1 almost everywhere in Q. If in addition Assumption 
\3.1\ holds we call the triple (v, ip, p) with the properties 

V G BC^{i- Lim n L\o, T- H^inf) , 

VeBC^{L;H\n))nLliJI;H'{n)), {^) e L\^{L- L\n)) , 
p£ Lli^^{I;H\n)) with Vpe L'^{0,T;L\n)) 
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a weak solution of ()3.ip - (|3.7p . if the following conditions are satisfied. 

- (/)v, dtilj)Q^ + (div(/)v v), + (2r?(v9)I)v, L>^)q^ - ((v ® J), 

= (/iV99,^)Q^ (3.10) 

/or a// e X (0, T))]'^ with divf/' = 0, 

- dtC)Q^ + (v • V^, Oq^ = - (m((^)V/i, VC)q^ (3.11) 
for allC(^C^mT);CH'^)), 

fi = ^''(93) + a'{ip)^- — ^ div {a{ip)Vip) almost everywhere in Qt and (3-12) 

(v,V^)lt=o = (vo,V5o) • (3.13) 

Moreover, 

Etot{^{t),v{t))+ f 2r]iip)\Dv\'^dix,T)+ f m((^)|V/i|2 d(x, r) 

<Etot{^{s)Ms)) (3.14) 

for all t G [s,oo) and almost all s £ [0, 00) has to hold (including s = 0). The total energy Efot 
is the sum of the kinetic and the free energy, see p.l6p . 

3.2 Existence theorem 

Our main result of this work is the fohowing existence theorem for weak solutions. 

Theorem 3.4. Let Assumption \3.1\ hold, vq G L^(il), ipQ G H^{Q) with \ipo\ < 1 almost 
everywhere and -f^ipodx £ (—1,1). Then there exists a weak solution (v, 99,^) of (|3.ip - (|3.7p in 
the sense of Definition [3!" 



In order to prove the theorem, we reformulate line (j3.12p to an equivalent equation. Therefore 
we introduce the function A{s) ■= a{T) dr. Then A'{s) = a{s) and 



-^/ai^AA{^) = a'{^) - div (a((^) Vip) 



2 



resulting from a straightforward calculation. 
With this notation (|3.12p reduces to 



/i = ^''(v9)- V^A^(V?). (3.15) 
We also rewrite the free energy with the help of A to 
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The kinetic energy is given by E-^^^{ip,v) = JqP-^ dx and the total energy as the sum of the 
kinetic and free energy 

Etoti^p, v) = E^^i^p, v) + Ei,M = X + X {^^"^^ ^ '^^2"^^'' ) (3-16) 

The next step is to rewrite (j3.15p with the help of a subdifferential. To this end, we set [a, b] := 
A([—l,l]) and define a reparametrized potential ^ through 

^ ' \ +00 else. ^ ^ 

This 4* fulfills analogous assumptions as that is ^' G C^{[a,b]) n C^((a,6)) and 

lim $'(r) = -00 , lim ^'(r) = +00 , $"(r) > -k, for aU r G (a, b) (3.18) 
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for some /5 G M. Here Assumption 13.11 (iv) is needed. We define ^o('') •= "^(r) + to get a 
convex function ^o- In particular it holds ^'{r) = ^'o(r) — /tr and ^''(j4(s)) = ^Q(yl(s)) — kA(s) 
for s G [—1,1]. Furthermore we observe ^{A{s)) = ^{s) for s G [—1,1] and therefore ^'{s) = 
^'{A{s)) • y/a{s), in particular (a(s))"2 . = 

This notation leads to a reformulation of ()3.15p to 

{a{^)r^ /i + kA{p) = %{A{ip)) - AA{^) . (3.19) 

Now we use a result from Abels and Wilke |AW07j for the energy E : L'^{n) with domain 
domi? = {ti G H^{Q,) \ a < u < b a.e.} given by 

Efu) = / ^"^^^'^^'^^^ + "^^2^°^'')'^^ foruGdom^, .^20) 
^ ' \ +00 else. ^ ■ ^ 

From Theorem 3.12.8 in |Abe07] . which is a variant of Theorem 4.3 in |AW07] . the domain of 
definition of the subgradient dE is given by 

V{dE) = {ue H^{n) I $'o(n) G L'^{n) , ^o{u)\Vu\^ G L^{n) , a„n|on = 0} 
and for u G T>{dE) it holds that dE{u) = — An + '^^{u). Furthermore there holds the estimate 

+ + X ^Hx))\Vu{x)\^ dx<C {\\dE{u)\\l, + + l) . (3.21) 

Analogously we define the energy E : L^{n) ^ R with domain domE^ = {p G H^{n) \ - 1 < 
if < 1 a.e.} given by 



E{p) 



( lf^\Vip\^dx + f^^o{^p)dx for^jGdomS, 
\ +00 else . 
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Here it holds that V{dE) = {if £ H'^{n) \ %{lp) G L^{n) , V(^p G L^{n) , = 0}. 

One can show that that A{ip) G T>{dE) if and only if G T>{dE) and therefore we get for 
u = A{ip) with if G V{dE) the identity 

dE{A{^)) = -^A{^) + ^'o(^M) . (3.23) 

This leads finally to the reformulation of (|3.12p as 

{a{^))-'^ li + KA{^) = dE{A{^)). (3.24) 

The fact that the right hand side equals the subgradient of E at the point A{ip) will be of 
importance in the following analysis. 

With the above notation we can rewrite the free energy as 

Eiree{v)= [ ^(v^) + / rfx = / $(A(^))dx+ / ^-^^^ dx 

= I ^o{A{^))dx+ [ ^^dx- / ~^{A{^)fdx 



We summarize the reformulation of problem ()3.ip - ()3.7p : 

dt{pw) + div(pv (g) v) - div(2r/((^)Dv) + + div(v (g) J) = /uVv? in Q, (3.25) 

divv = in Q, (3.26) 

(9tV3 + V • Vv? = div {m{ip)Vti) in Q, (3.27) 

(a(v9))~^ /i + K^(v9) = ^'o{A{ip)) - AA{^) in Q, (3.28) 

^r\^Q = on 5, (3.29) 

dn^lan = dniAdn = on S", (3.30) 

(v,99) |t=o = (vo,<^o) in ^, (3.31) 

where J = — ^^2^^ m((^)V/i. 



4 Implicit time discretization 

4.1 Definition of the time-discrete problem 

In order to prove Theorem 13. 4| we use an implicit time discretization. To this end, let h = for 
iV G N and G (Jl), G //^(Sl) with ^''((^fc) G L2(J7) and = \{pi + P2) + \{p2 - Pi)^k 
be given. We construct (v, 99,/i) = {vk+i,ipk+i, Pk+i) as solution of the following non-linear 
system, where 

J = Jfc+i = -£^m((^fc)V/ife+i = -£2^m(99fc)V^ . 
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Find (v,^,/i) with V G H^{nfnLl{n), (f G V{dE) and ^l G HHQ) = {n G | dnu\g^ = 

on 517}, such that 

pv-^PfcVfc ^^^ + (div(/>fcv » v), + (27?((/^fc)Z)v, Z)V)n + (div(v » J), V') ^ 

= (/iVv9fc,t/»)^ (4.1) 

for all -0 G Co";j,(r?), 

^ + V • = div {rn{Lpk)V ^) almost everywhere in and (4-2) 

~ '^^ + K ^'^'^^ + ^(V^fc) ^ _AA((^) + ^ol^lv')) almost everywhere in Q.. (4.3) 



^((/.)-^(^fc)' 

Remark 4.1. (i) Multiplying identity ()4.2p mt/i — P'^^p^ leads to 

P- Pk 



h 



V • Vpfc = div J , 



where we used the definition of J = 3k+i o^nd the linear dependence p{(p) = 2(pi + P2) + 
^(P2 — Pi)v^- Using div(v (8) J) = (div J)v + • v leads to an equivalent version of 
()4.ip given by 



h 



P''^^ ^) +(div(pfcV»v),0)^ + (2?7(v?fc)Z)v,Z)^)f, 
n 



+ ( (div J - _ V . Vpfe) I , ^ + ( ( J • V) V, V) ^ = il^ViPk , (4.4) 

for all xjj G C^^{Vl). We will use this equivalent version in the following especially in the 
a-priori estimate for solutions of the time- discrete problem. 

[ii) Integrating equation ()4.2p with respect to the spatial variable, using div v = and the 
boundary conditions, we obtain J^ipdx = ipk dx, which means that ipk dx = tpo dx 
is constant. 

Lemma 4.2. Assume that ip G T>{dE) and /i G H^{Q,) solve ()4.3p with given ipk G -ff^(il) with 
\fk\ ^ 1 in Q such that 

1^ y (/? = y Lfkdx £ (-1, 1) . 

Then there is a constant C = C{J^ (pk) > 0, such that 



l^o(^M)llL2(f^) + 



/i dx 



n 



< Ci\\Vfih2 + \\Vp\\l, + WVpkWh + 1) and 
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Proof. First we note that F{ip,ipk) := fulfills < c < 99^) < C for some c, C 

independent of A; due to 



A'{Tip + {l-T)ipk)dT = / \J a{TLp + (1 - T)ipk)dT 



and < mo < a, a G C"'^(IR) and (p{x) G [—1,1] almost everywhere. Testing ()4.3p with ^ 
ipk){<P> — where ^ = -py /^^ is the mean value of (p we get 



-AA(V9)F(v9,(/Pfc)((/?-^)dx+ / $o(A((^))F(v7,(^fc)((^-^)(i2;. 



(4.5) 



With HQ = fi — fj, we observe that fi{ip — dx = fiQ^ dx. Furthermore, due to 



= (^a-^T^ + il-r)ipk)a{Tip + {l-T)^k){rVip + {l-T)V^k))dr, (4.6) 
the boundedness of a > niQ > and a' on [—1, 1] and (r(/? + (1 — T)ipk) G [—1, 1], we get 

<Ci\Vip\ + \V^k\), 



^/ Ai^)-Ai^k) \ 



which leads to 



-AA{ip) ■ F{ip, ipk){(p - ^p) dx 



Jn ' — V — ' VV 'P-^k J 



{ip — ip) \ dx 



I a{Lp)'SIip 



<C [ \Vip\^ dx + C [ \Vip\ (|V(^| + \Vipk\) dx 
Jn Jn 

<C{\\V^\\l, + \\V^k\\h) . 



To estimate the last integral in (|4.5p we calculate at first '$q{A{ip)) = ^0(92) + KA{ip) — Kip 
and use the fact that Tp E (—1 + e, 1 — e) for some e > 0. In addition with the assumption 
lim^_j.±i ^'o(</') = ±00 one can then show the inequality '$Q{ip){ip — Tp) > C£|^o((/?)| — in three 
steps in the intervals [— 1, — 1 + |], [— 1 + |, 1 — |] and [1 — |, 1] successively. Altogether this 
leads to an estimate of the following form: 



%{A{cp))F{ip, ipk){ip-^)dx>C \%{ip)\ dx-Ci. 
n Jn 
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Using the last inequalities and ()4.5p we get 



\%{^)\dx<C\\^iQ\\L2y\\L2+C / -\A{^) + A{^k) 

Jq ^ ' 



<ciM+\vk\) 



A{^) - A{ipk) 



\ip — ip\dx 



+ Ci\\V^\\l2 + \\V^k\\h) + C, 
<Ci\\^io\\L2 + \\V^\\l, + \\V^k\\h + l) 
<C{\\Vf,h2 + \\Vip\\l, + \\V^k\\l2 + l). 

where we have use |^((/')| < C\ip\ and the bounds \ip\, \ifk\ < 1- Now we test ()4.3p with 
to get 



A{ip)-A{ipk) 



Hdx= VA{ip) ■ V 



n 



f-^k 



which together with the previous estimates leads to 



fi dx 



<C(||VMk2 + ||V(^||i. + ||Vv9fc||i. + l). 



Finally, the estimates of the subdifferential dE{A{Lp)) = -AA{lp) + '^'q{A{lp)) and of ^'0(^(9?)) = 
^'q ((/?), each in the L^(J7)-norm, follow directly from ()4.3p and inequality (|3.2ip . □ 

4.2 Existence of a solution of the time-discrete problem 

Lemma 4.3. Let Vk € L^(r2), (pk S H'^{^) and pk = ^{pi + P2) + \{P2 — Pi)^k be given. Then 
there are some iv,ip,p) G {H^{n)'^ n Ll{n)) x V{dE) x H^{n) solving (glD-gSl), which satisfy 
in addition the discrete energy estimate 



Etotiv, v) + / Pk ^ dx + 



\VA{^) - VA{^k) 



dx 



+ h 2r]{ipk)\Dv\'^ dx + h / m{ipk)\Vp\'^ dx < Efotiv'k^^k) ■ 



(4.7) 



Proof. First we show the a-priori estimate (fiTT]) for any (v, ip, p) G [Hl{rtY ^ Ll{Vi)) x V{dE) : 
Hl{^) solving (liTTD -M. 

In order to test ()4.ip with i/) = v we need some preparations. First we observe 



(divJ)-+ (J-V)v) •vdx= / div J 



dx = id 
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and then we calculate 

" dx 



j (div(pfcv v) - {Vpk • v)^) • V ( 

= j (^div(/)fc V (g) v) - div(/9fc v) • v dx 

div Pk^—^ dx = , 



where we used integration by parts and divv = 0. Furthermore from the simple algebraic 
equation 

a.(a-b) = ^-^ + ^^^ fora,bGM'^ 

we get 

^ (pv - pfeVfc) • V = i(p - Pfc) |v|^ + ipfc (V - Vfc) • V 

n n n 

_ M__i2 , 1 _ fH^ Wk\^\ , 1 _ |v-Vfcp 



= ^(P-P.)|v| ^j+X^'^ . 

Testing ()4.ip or equivalently ()4.4p with tp = v and using the above identities we obtain 

0= / /" /" 2r?(v?fc)|I)v|2dx- /" ^ (V(pfc • v) dx . (4.8) 

Moreover, choosing p in ()4.2p as a test function, we get 

= / p.dx+ / (v- V(/?fc)^dx+ / m(v?fc)|Vp|^(ix. (4.9) 

Finally we test (|13D with - M'^k)) to get 

=\ ^ VA(^) • V{A{^) - A{^,)) dx + J^ q.'^(A{^))^^^±zAk^ dx 



,^dx- [7i^^^)^-^dx. (4.10) 
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Summing the identities ()4.8p - ()4.10p leads to 

Jo. n Jn ^iT' 

+ i ^ VA(^) • V{A{ip) - A{^k)) dx 

> I — — f^^^^^ dx+ I pk ^ dx+ I 2-r]{ipk)\Dw\^ dx + [ m{ipk)\V pf dx 

1 f |V^(y) - VA{nr I 1 f /EMv)l_\YMVi)l)j:,^ 
h Jq 2 h Jq \ 2 2 J 

where we have used 

%{A{^)) {A{ip) - A{ipk)) > ^o{A{v)) - ^o{A{^k)) and 

v^M . v(^(.) - AM) ^ M _ ^iM^ + 1^-^'^' -/-^'^^"^ . 

This leads to the claimed discrete energy estimate given by 

Etot{f,^)+ / Pk dx+ / dx 



Jn ^ Jn ^ 

+ h 2ri{(pk)\Dv\'^ dx + h / m{(pk)\^ dx < Etot{^k,^k) ■ 
Jn Jn 

In order to show existence of weak solutions we want to use the Leray-Schauder principle. With 
the abbreviation H^{^) = {u G H'^{^) \ dnu\Qn = on (90} we define operators Ck,J^k ■ X — )■ 
y, where 



X = [Hli^Y CMli^)) xV{dE)xHl{n), 

Y = (^H^{nf n Ll{n)y X L^{n) x L^{n) . 

For w = (v, if, n) G X we set 

A(w) = ( -div(m(99fc)V^) + j^pdx 
A{^)+dE{A{^)) 



where 



(Lfc(v),^)= / 2r]{ipk)Dv.D'ipdx for iP e H^{n)'' n lKQ) 
Jn 
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and the second and third hne are regarded pointwise. Note that due to (p G T>{dE) it holds 
A{(p) £ and therefore the last line in £fc(w) lies in L^(il). Furthermore for w = (v, ip, fj.) G 

X we define 



/_pv^ - div(p,.v 0^r)+ fiV^k - (div J - £z£^ - v • Vp^) f - ( J • v) v\ 



h 



\ + + 2 / 

Then w = (v, (p, fi) G X is a weak solution of the time discrete problem (|4.ip - ()4.3p if and only if 

£fc(w)- J-fc(w) = 0. 

Note that we used the equivalent version ()4.4p instead of ()4.ip . that is we substituted div(v(g) J) 
accordingly. 

From standard theory of partial differential equations we get the invertibility of 



Now we consider for given / G L^(0) the elliptic boundary value problem 

-dw{m{(pk)V fi) + / fidx = f inn, dnfi\QQ = 0. (4.11) 
Jn 

With the Lemma of Lax-Milgram we get the unique existence of a weak solution /i G H^{Cl) and 
we want to show with a bootstrapping argument that this solution fulfills fj, G H^{Q,). At first, 
we observe that fi is also a weak solution to the following problem: 

An = - {m{ipk)y^ (v{m{ipk))-Vfi + j^fidx- inn, dnfi\Qn = 0. (4.12) 

Since ip^ G H'^{n) with < 1 almost everywhere in n, we conclude V(m((^fc)) G L^(il) and 
for the product of an L^- with an L^-function we get V (m{ipi^)) ■ Vfj, G (n). Using the bound 

~ 3 

of m{ipk) from below by a positive constant, we get for some / G L'2{n) that 

A/i = /GLi(Jl). 

Elliptic regularity theory now gives /i G W'i{n). In particular this means G 1^3(^2), which 

2 ^ 2 

embeds into L'^(O), so that we conclude for the product of an L^- with an L'^-function V {m{ipk)) ■ 

V/i G L2(0). This shows 

A^l = jeL\n), 
and therefore finally /i G H'^{n) together with an estimate 

II//IU2 <C(||/.||^i + 11/11^2) . (4.13) 
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Since an i?^-solution of ()4.12p also solves ()4.1ip . we obtain that the solution of ()4.1ip lies in 

The operator dE is maximal monotone and hence 

I + dE : V{dE) L^(Sl) 

is invertible. In addition we have continuity of the inverse considered as a mapping L^(f2) — >■ 
H'^~'^{Q) for arbitrary < s < j. This is shown in the proof of Proposition 7.5.5 of Abels |Abe07] 
for a similar operator. For the convenience of the reader we give the details. Let /;—)■/ in 
L^(ri) for / — )• c« with fi = ui + dE{ui) and f = u + dE{u) be given. Then we have — )■ n in 
H'^{n) since 



\ui 



u 



1^2 + llVuz - Vu|||^2 < \\ui-u\\i2 + [dE{ui)-dE{u),ui 



< \\ui + dE{ui) 

<^ll/^-/lli2 + ^lh 



{u + dE{u))\\L2 \\ui -u\ 



L2 

L2 



Moreover, because of the estimate ()3.2ip . ||n;||j|^2 is bounded and together with compactness we 
get the convergence tt; — )■ n in for arbitrary < s < \. 

Since A is invertible and we have the equivalence A{^p) G T){dE) if and only if 99 € 'D{dE), 
we also see that 



A{.) + dE{A{.)) : V{dE) L^i^) 

is invertible with continuous inverse considered from L'^{Q.) — )• H'^~^{Q) as above. Altogether we 
obtain that : X ^ Y is invertible with inverse C^^ : Y ^ X. Note that X is not a Banach 
space since T>{dE) includes inequality constraints. 

To get a continuous and even compact operator, we introduce for < s < ^ the following 
Banach spaces 

X := [H^in)'' n l2 (17)) X H^-^{n) x Hl{n) , 

Y := Li(J])^ X Wl{n) X H^{n). 

2 

Then we obtain continuity of : Y ^ X from standard theory and with the above note 
concerning the continuity of the third line in C'f^^ . 

Finally due to Y Y the restriction -.Y^Xisa compact operator. 

The next step is to observe that Tk ■ X ^ Y is continuous and maps bounded sets into 
bounded sets. More precisely we obtain the following estimates (note that (pk and therefore pk 
lie in H^{n)): 



||/0V||^3 <C||v||Hl(||(^||i2+l), 
\\fJ.V<Pk\\^3 < Ck\\fJ.\\L2 , 

||(J-V)v||^3 <C||v||hi||/.||^2, 



||div(/)fcV(g) v)||^3 <Cfc||v|||i, 
||(div J)v||^3 < Cfc||v||^i||^||^2 
l|v • V(^fc|liyi < Ck\\v\\H^ , 



\\Ai^)\\Hi<CM 



< Cfc||^||//2(||y?||^l + 1) 
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In detail, these eight bounds are derived as follows: 

(i) For pv, the estimate is clear due to v G ^ L^. The appears, since p depends 

affine linear on ip. 

3 

(a) For the estimate of div(/3fcV (S) v) in L2 , we have to estimate on the one hand terms of the 
form pi^diViVj in L2 , which are a product of functions in L°°, and and therefore are 

3 3 

bounded in L2 . On the other hand we need terms of the form dipkVi^j in L2, which are 
a product of three functions in L" and therefore are bounded in particular also in L2. 

{iii) The estimate for fiVipk in L2 follows immediately from p £ and diip^ £ L^. 

(iv) For the estimate of (divJ)v in , where J = — m{ipjS)V p, we have to estimate on 
the one hand terms of the form ■m'{ipk)di(pkdjpvi, which are a product of functions in L°°, 
L^, and and therefore are bounded in particular in L2. On the other hand we have 
to estimate terms of the form m{ipk)didjpvi, which are a product of functions in 
and L" and therefore are bounded in L 2 . 

(v) The estimate of (J • V)v in L2 follows since the terms m{ipk)dipdjVi are a product of 
functions in L°°, and and are therefore bounded in L2. 

(vi) For the derivatives of v • V^pk, we observe on the one hand that 9jV; is bounded in and 
dj^pk is bounded in L^. On the other hand Vj is bounded in and didjip^ is bounded in 
L^, so that in both cases the product is bounded in L2. 

(vii) The estimate of A(ip) follows from |^((^)| < K\ip\ and di{A{ip)) = A'{<p)diip almost every- 
where, where A'[(p) = a{'p) < K. 



{via) To estimate 



p, we use first 



-1 



A{^) - A{ipk) 



\l aijip + (1 - T)<Pk) dr 



and then from ()4.6p in the proof of Lemma 



>c>0 

the result 
<C(|V(^| + 1) 



1 

< - 

c 



In particular this gives a bound on 



V 



A{ip)-A{ipk) 



V-^k \ ^ f A{ip) - A{ipk] 
A{v)-A{ipk)) V H^-^k 



Together with p € H'^{Q) ^ C^{Q) we can bound V 



A{p)-A{ipk) 



p in and of course 



also 



Vp and 



p each in 
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We want to use the Leray-Schauder principle on Y and rewrite therefore the identity £fe(w) — 
J),,(w) = for a solution w G X of (jil^ - ijO]) into 

f- J-feo£^i(f) = for f = £fc(w). (4.14) 

We set /Cfe := J^fco£^^ -.Y and note that /C^ is a compact operator because is compact 
and J-fc is continuous. Equation ()4.14p is then equivalent to finding a fixed point of /Cfe, that is 

f-/Cfc(f) = ^ f = /Cfc(f). 

To deduce the existence of such a fix point with the help of the Leray-Schauder principle (see 
for example Zeidler |Zei92| ). we have to show that 

> such that, if f G y and < A < 1 fulfills f = A/C^f , then ||f||y <R. (4.15) 

So let f G y and < A < 1 with f = A/Cfcf. With w = ^(f) we have 

f = A/Cfc (f ) ^ Ck (w) - A J-fc (w) = , 

which is equivalent to the following weak formulation 

/ 2r){(pk)Dv : D'il)dx + X [ ~ P^^^ -ipdx + X [ diY{pkV (g> ■ tp dx 
Jn Jn " Jn 

+ X j ^divJ - ^~ - V • Vpk^ ^-ipdx + X j (j -V^v ■ tpdx = X j fi\7(pk ■ ip dx 

(4.16) 

for all tp G H^{nf H LKQ) and 

A^— — ^ + Av • V(/3jfc — A / /i = div(?Ti((/7,fc)V/x) — / fidx, (4-17) 
h Jn Jn 

Ai^) + dEiAi^)) = XAi^) + ^J^^^^^f + A7^^iM±^ . (4.18) 

First we derive an estimate for w = (v, ip, p) in the norm of X and additionally an L^-estimate 
of dE{A{ip)), then we conclude the desired estimate for f due to the boundedness of Fk- 

Analogously as in the derivation of the energy estimate ()4.7p we set ■?/> = v in ()4.16p . test 
()4.17p with /i and (|4.18p with j^{A{ip) — A{ipk)) to get with similar calculations (we omit the 
integration element dx for reasons of shortness): 

+ / m(v9fc)|VH' + (l-A)i / A{^){A{^)-A{ipk))+'^ j VA{ip)-{VA{^)-VA{ipk)) 
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2 2 

|2 



+ / m{^k)\Vfi\' + {l-X) . , 
In 'T' JQ \ 



2 /l/,^ \2 



1 



|Vyl((/PA 



This leads to the fohowing estimate 

h I 2r,{^k)\D^? + h [ m((^fc)|V^|2 + i [ \VA{^)\^ + [ ^ {A{ip)) + {I - \) 
JQ Jn ^ Jn Jn 

< [ ^ + 1 [ Ai^,f + l [ |VA(^.)P + 



/X dx 



n Jn 



where we omitted the (nonnegative) terms A dx, A Pk^-^^—j^ and (1 — A) "^^^ 
on the left side, since due to the factor A resp. (1 — A) they will not give a contribution to some 
estimate of ||w|| ^ independent of A. Note that due to w = (v, (p, 12) = £,^^{f) G X, it holds that 
if € T)[dE) and therefore (p G [—1,1] almost everywhere, which implies in particular /9 > 0. From 
this fact we also get boundedness of the term ^[A{ip)) dx, which can be estimated therefore 



on the right side. We also used the simple estimate — AJ^k 2'° dx < X \k\ ^ dx and in 
addition we estimated every A resp. (1 — A) on the right side against 1. 
Using |Vj4((^)P = a((/3)|V93p, we summarize the previous estimate to 



(1-A) 



P 



dx^ +h f 2rj{ipk)\Dv\'^ dx + h f m((/?fc)|V^p 
/ Jn Jn 



dx ^ — / a{ip)\V(p\'^ dx < Ck ■ 
2 Jn 

(4.19) 



With the help of < C due to (p G V{dE), Korn's inequality for v G Hl{^)'^ r\ L'l{^) and 

the fact that r/, m and a are bounded from below by a positive constant, we get the estimate 



In 



/i dx 



+ l|v 



(4.20) 



To get an estimate of the L^-norm of the chemical potential ^, we differ two cases. If A G 1], 
then we proceed with the simple estimate ^| j^ndx\ < A| f^fj.dx\ and get as in the proof of 
Lemma 14.21 together with ()4.20p from equation ()4.18p the inequality 



/.i dx 



<Ck. 



For A G [0, ^) we use (|4.20p directly to get also here \J^^dx\ < Ck- With this inequality for the 
mean value of the chemical potential we can improve estimate ()4.20p to 



(4.21) 
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Together with ()4.13p we also get an estimate of the //^-norm of the chemical potential /i given 

by 



+ Il/^ll/f2(n) + llv'll/fiCQ) < Ck ■ (4.22) 

Using identity (|4.18|) pointwise we also get ||(?£'(A((/7))||^2(q) < Ck- In summary this leads to 

||w||^ + \\dEiAiipmL2^n) = ||(v,<^,m)||^ + \\dE{A{^))\\L2^^^ < Ck. 

To get finally an estimate of f = >Ca:(w) in Y we use that f — A7"fc£^^(f) = implies f = AJ^fc(w) 
and the fact that Fk '■ X ^ Y maps bounded sets into bounded sets, which holds due to the 
above estimates for Fk- This gives 

||f||y = ||AJ-fc(w)||^ < Cfc(||w||^ + 1) < Cfc , 

which was the remaining part to apply the Leray-Schauder principle as described above. There- 
fore we proved the existence of a weak solution of the time discrete problem (j4.1|) - (|4.3p . which 
additionally fulfills the discrete energy estimate ()4.7p . □ 



5 Proof of the Main Result 
5.1 Compactness in time 

To complete the proof of Theorem 13.41 we have to pass to the limit /i — t- resp. — )• cxd in 
our approximate solution. Therefore let S N be given and let {^k+i-,^k+i-, ^^k+l) be chosen 
successively as a solution of ()4.ip - ()4.3p with h = and (vq, ^Pq) as initial value. Here we have 
to approximate the initial value lpq S H^{Q) by i?^-functions Lp^ in order to apply Lemma 14.31 
to the first step. This is done with standard partial diff'erential equation theory. For example 
one could choose u as the solution of the following heat equation 



dtu — Au = 





in 


n X (o,r) , 




u = 


fo 


on 


nx{t = o}, 









on 


dn X (0, T) 






\^^\ 


< 1 


in J7 and - 


■> ifo in H^{n) 



and set ipf^ := u\,_ i to 

Then we define f^{t) on [—h,oo) through f'^{t) = fk for t G [{k — l)h,kh), where A; G No and 
/ € {v, If, /i}. Also set /o^ = ^{pi + P2) + |(P2 — Pi)p^ ■ In particular it holds f^{{k — l)h) = fk, 
f^{kh) = fk+i and f^{t) = fk+i for t G [kh, {k + l)h). Additionally we define 

(A+/) (t) := fit + h)- fit) , (A^/) it) := fit) - fit - h) , 
d^Jit) ■■= I i^tf) it) > d-Jit) := I (A-J) it) , 
h:= irU) it)=fit-h). 
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Then for arbitrary ijj G (C^(r2 x (0, oo)))'^ with divi/> = we choose ij) := jj^^^^ ijj dt as test 
function in (|4.ip and sum over G No to get 



poo p poo p poo p 

/ d^^{p^w^)-xpdxdt+ div{p^v^ (g)v^) ■tpdxdt+ 27j{^p^)Dv^ : Dtp dx dt 

Jo Jn ' Jo Jn Jo Jn 



(S)J^) : Dipdxdt= I p^Vif^ -ijjdxdt (5.1) 
JQ. ' Jo Jn 

for all G (C^(r2 x (0, oo)))'^ with divtp = 0. The first term can be rewritten due to 

poo p poo p 

/ / (P^^^ )-ipdxdt = - / ) • dx dt . 

Jo Jn ' Jo Jn ' 

Analogously we get 

poo p poo p poo p 

/ / d-,^ip^ C dx dt + / / v^^f^ ■VCdxdt= / / m{^^)V p^ ■ V C dx dt (5.2) 
Jo Jn ' Jo Jn Jo Jn 

for all C e C^{{0, oo); C^{n)) and 

- N ~ 

■p'' + UMv>'')+A{^^)) = -AA{^^) + ^',{Ai^^)) (5.3) 



A-^((/,^) 



holds pointwise in $7 x (0, oo) almost everywhere. 

Let [t] be the piecewise linear interpolant of -E'tot(9'fc; v^,) at t^ = kh given by 

(k + l)h — t t — kh 

E^{t) = ^ h Etot{^k,^k) + —^Etot{Vk+i,^k+i) for te[kh,{k + l)h). 

Also define for all t £ (t^jifc+i), k £ Nq 

D^{t):= [ 2r?((^fc)|Dvfc+i|2dx+ [ m{^k)\^ fik+i\^ dx . 



Jn Jn 
Then the discrete energy estimate ()4.7p implies 

-—E^{t) = -^tot(yfc, Vfc) - -Etot(v?fc+i, Vfc+i) > ^^Af^^) (-5 4) 

for all t G (tfc)^fc+i)) G Nq. Multiplying this inequality by r G ^^-^'-'^(0, 00) with r > 0, 
integrating and using integration by parts gives 

poo poo 

i^tot(<^^,vo)T(0)+ / E''{t)T'{t)dt> D''{t)r{t)dt. (5.5) 
Jo Jo 

Integrating (|5.4p gives 

i?tot(v^^(t),v^(t))+ r/ (2,7(^;y)|Z)v^|2 + m((^j^)|V/.^|2)dxdr 
Js Jn 

<E,^,{^^{s\^^{s)) (5.6) 
22 



for all < s < t < oo with s,t £ HNq. 

Together with Lemma [^^2] and the fact that Stotlv'o^; ^o) is bounded this leads to the fol- 
lowing bounds: 



is bounded in L'^{0,oo; H^{nf) and in L°°{0,oo; L'^iQ^) 
V/i^ is bounded in L'^{0,oo; L^{n f ) , 



N 



is bounded in L'^{0,oo; H^{n)) and 



(5.7) 



for a monotone function C : 
get 



So \In ^^1 dt < C{T) for ah < T < oo 

— >■ M+. Using these bounds, we can pass to a subsequence to 



V 



N 



N 



V in L'^{0,oo;H^{n)'^), 

in L'^{0,oo;L\nf)^ (^L\0, oo; L^flf) 

^* if in L°°{Q,oo;H^{VL)) ^ {L^{id,oo\H^{^)))' 
/i in L'^{0,T;H^{n)) for ah < T < oo , 
in L'^{0,oo;L'^{nf) . 



Here and in the following all limits are meant to be for suitable subsequences — )■ oo (resp. 
/ifc — )• 0) for /c — oo, unless otherwise stated. 

Now let ip^ be the piecewise linear interpolant of ip^{tk), where tk = kh, h G No, i.e. 



'N 



^X[o,/i] *t J where the convolution is only taken with respect to the time variable t. 



Then it holds that dtf 



■N 



\h'P and 



N w ^ 1, 1 1 n ~l\ 1 1 



AT I 



(5.8) 



From line ( 15.211 we get that dtif^ G L2(0, oo; ^-^(f])) is bounded, since v^(/j^ and V/i^ are 
both bounded in L^(0, oo; L^(i7)'^). Together with boundedness of (p^ in L°°(0, oo; //""^(il)), 
which follows from the boundedness of ip^ in L~(0,oo;i?^(il)), we get with the help of the 
lemma of Aubin-Lions (j2.ip the strong convergence 



ip in L'^{0,T;L'^{n)) 



for all < T < oo for some p £ L°°(0, oo; L^(il)). In particular it holds for a subsequence that 
ip pointwise almost everywhere in (0,oo) x il. Additionally by the above estimate ()5.8p 



■N 



"P 

we have 



■AT 



'P 



N 



in L'^{0,oo-H-^{n)) 



which gives p = p. Furthermore, since p^ G H^^^^{[0,cc); H ^(f^)) n L^j^^([0, oo); i/^($7)) ^ 
BUC{[0,oo);L'^{n)) and p^ E L°°{0,oo; H^{n)) are bounded, it follows from Lemma O that 

ip£ BCU[o,^y,H\n)). 
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Finally it follows from the bound of in H\0,T; H~Hn)) and in L°°(0,T; for 
< T < oo together with the fact that 99 in L^{0,T; L^{n)) that ^^(0) ip{0) in 

L'^{Q). But the left side equals 99^, which converges to ipo in L^(J7) so that we finally conclude 
(f{0) = (fo. 

Analogue observations can be done for since it depends affine linear on 99^. 

To show the convergence of ()5.3p . we observe that the right side is given by dE{A{ip^)) and 
this is bounded due to Lemma 14.21 in L^(0, T; L^(il)) for < T < 00. Moreover, the left side 
converges weakly in Lp'{Q,T] L'^{Vt)) to 

/ ■■= a{ipy^^ p, + RA[(p) , 

which means that dE{A{ip^)) f weakly in L^(0, T; L'^{Q)). If we now show that 

limsup(5^(A(^^)),yl(<^^)) < {f,A{^)) , (5.9) 

we can use the fact that dE is a maximal monotone operator and apply Prop. IV. 1.6 in Showalter 
|Sho97j to conclude that dE{A{Lp)) = f, which would finally lead to (j3.12p . 

But the result (|5.9p does hold even with equality, since with the natural definition of as 
the left side of ()5.3p . we have 

= (/^,A(^^)> {f,A{^)) 

due to the strong convergence of ip^ in L^(0, T; and therefore also of A{ip^). 

Now we use the estimate ()3.2ip for A{ip^), which together with Lemma 14.21 gives bound- 
edness of A{ip^) in L^(0, T; Then the same holds true for ip^ and due to the strong 

convergence of (p^ in L^(0, T; we conclude with an interpolation argument even the 

strong convergence 

ip^^if in L'^{0,T;H^{n)). (5.10) 

The next step is to show strong convergence — )■ v in L^(0, T; L^(ri)°') for all < T < 00 to 
conclude a convergence pointwise almost everywhere. As above let pv^ be the piecewise linear 
interpolant of (p^v^) (t/.), where tk = kh, h G Nq. Then it holds that dt = d^j^ [p^v^). 

With the help of the projection onto L'^{^) we get from line (|5.ip that dt ^Po-(pv^)^ is 

bounded in Lr (0, T; W^:^ (17)) for < T < 00, since due to some known interpolation inequalities 
we have the following bounds 

p^v^(^v^ is bounded in L^{0,T; d {Q)) , 

Dv^ is bounded in L^{0,T; L^{n)) , 

v^®V/i^ is bounded in lI {0,T; lI {Q)) , 

p^Vph is bounded in l2(0, T; Li(0)) . 

In detail, these bounds are derived as follows, where for abbreviation "E" always means bounded 
independent of N in the corresponding space and we omit the time and space variables. 
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{i) Due to G L°°{L'^) n L'^{L^), we get for products vf vj^ G ^^(Lt) and from £ 
L°°(1^t) this also holds for p^v^ (S) v^. 

(iz) From v G L'^{H^) we immediately get Dv^ G L^(L^). 

(m) Due to G L°°(L2) n L'^{L^) and V^^ G ^^(L^) we get vfdjfi^ G ^^(Li) n L\lI). 
Now, if / C M is an interval, {Xo,Xi) is an interpolation couple of Banach spaces, and 
X = {Xq, Xi)g^q is the real interpolation space of type 9 G (0, 1) with I < q < oo, we have 
the embedding 

LP^{I;Xo)nLP'{I;Xi) ^ LP{I;X), where - = + — , (5.11) 

P Po Pi 

where 1 < po,pi < oo. In addition, we use the description of real interpolation spaces with 
< ^ < 1 for Lebesgue-spaces with 1 < qo,qi < oo given trough 

(L<?o(J7),L'?i(S7))0„ = where - = i^ + — . (5.12) 

For Xq = L^{Q,) and Xi = L^{Q) we get with = j that <? = | and therefore X = 
{Xq, Xi)g^q = L3{Q). Furthermore with po = 2 and pi = 1 we get p = f , which leads to 
the bound vfdjp^ G L?(Lt). 

(iv) From /i^ G L^{L^) and V(y9^ G L°^(L2) we get finally p^Vif^ G ^^(Li). 

This means in particular that all four terms are bounded in L7 (0, T; La (0)). Therefore we can 
allow in (|5T]) for test functions with ^, G (^L?(Lt)^ = L^{L^), i.e. ^ G -L^(W2). This 

implies the bound for dt (p<,(p^^)) in {L^{Wl))' = L^W^^). 

Additionally P<,(pV ) G L^{0,T; H^n^) is bounded and we can therefore conclude with 
the help of the Lemma of Aubin-Lions (j2.ip the strong convergence 

P,(^^)^w in L\0,T;L^{nf) 

for ah < T < oo for some w G L°°{0,oo; L'^{n)'^). 

Since the projection F„ : L^(0, T; L^(J7)'^) — ^ L^(0, T; L^(r2)) is weakly continuous, we con- 
clude from the weak convergence pv^ pv in L^(0,T; L^(ri)) that w = ¥fj{p'v). 

Now we derive the strong convergence — )■ v in L^(0, T; L^(il)'^) with the help of the 
following observations: 

r [ pV\'= f I P.(p^v^)-v^^ f [ P.(pv).v= H' pH\ 
Jo Jn Jo Jn Jo Jn Jo Jn 

where we used the strong convergence of Pct(/9^v^) in L^(0, T; L^(il)'^) and the weak convergence 
of in L2(0,r;L2(J7)rf). This gives (p^)5v^ ^ (p)5v in L'^{0,T; L^i^lf) and from above we 
know that 

p'^ —?■ p almost everywhere in (0, oo) x Q and \p^\ > O , 
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so that we can conclude 

v^ = (p^)-^ ((p^)iv^) in L\0,T;L\nr). 

This means in particular that — t- v pointwise almost everywhere in (0, oo) x $7 (for a sub- 
sequence). We note that similar ideas to prove strong convergence of have been used by 
P.-L. Lions |Lio96t Section 2.1]. 

Using these convergence results together with the fact that for all divergence free ■0 the 
following convergence holds 

Jn Jo Jn Jo Jn Jo Jn 

we can pass to the limit in the equations (jS.ip . (j5.2p to get (j3.10p . (j3.1ip . 

5.2 Initial data for the velocity v 

Lemma 5.1. Let v, v £ Ll{il.) D H^i^l.)''- and p £ L°°{il.) with p>c>0 such that 



/ pv -tpdx = / pv ■tjjdx for all ip G C^ai^) ■ 
Jn Jn 



Then it holds that v = v almost everywhere in Q. 

Proof. By approximation the identity also holds for := v — v E L'^{^) H Hq{Q,) and we get 
JqPoW ~ vpdx = 0, which due to the boundedness of po from below by a positive constant 
gives finally v = v in L^(il)'^. □ 

With the following arguments we show that the velocity v fulfills the initial condition (j3.13p . 
At first we derive weak continuity in time for the projection of pv onto the divergence-free vector 
fields. Therefore let := Fa-{pv^) and from the arguments above we deduce the boundedness 
of 

in Wl^^^J[0,<x^);W^\n)) ^ BUC{[0,^);W^\Q)) and 
in L°°{0,oo;L'^{nf). 

Due to these bounds, we concluded above that — ?• w in L^{0,T;L^{nf) and with the help 
of Lemma [2TT] we deduce for the limit that w € BCw{[0, oo); L'^ (0,)^) . Now let < T < oo and 
consider the auxiliary problem for a function p £ L'^{0,T; H^q^{Q)) given by 

-div(^V^,(t)) = div(^w(t)) in n, 

Vp{t) • n = on (90 . ^ ' 

The corresponding weak version reads as 

- I ^Vp(t)-VTdx= I ^w(t)-VTdx for all TeH\n) (5.13) 
Jn p{t) Jn p{t) 
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for almost all t in (0,T). By the Lemma of Lax-Milgram there exists a unique solution and it 
fulfills the estimate 



\\^p{t)\\LHn)<C\\^mLHn), (5.14) 

where C is independent of t. Here we used the fact that p G L°°{Qt) is bounded from below by 

a positive constant. 

Claim: Vp £ BCU[0,T]; L'^{nf). 

To see this claim, let tn, t G [0, T] for n G N with t„ — t- t be given. 

we conclude that in particular w(t„) is bounded in and from the last estimate (|5.14p also 
\7p{tn) is bounded in L^. Therefore we get the weak convergence Vp(tnJ ^ G L^{n)'^ for 
A; — )• oo at least for a subsequence, since {Vq \ q G H^qs^{U)} is a closed subspace in L^(ri)'^. But 
due to the unique solvability of (j5.14p we arrive at Vq = Vp{t). Since this argumentation holds 
for any weakly convergent subsequence, we even have Vp{tn) Vp(t). 

In addition we already know due to the weak continuity of Pq- and due to the weak conver- 
gence pv^ — )• pv in L^(0, T; L^(il)'^) for every < T < oo, that w = Po-(/9v). In particular this 
gives 

w(t) = Po- (p(t)v(t)) almost everywhere in (0,oo) . 

By characterization (|2.2p of the projection Po- onto L^(il), the last identity means 

w(t) = p{t)v{t) — Vp{t) almost everywhere in (0,c«) , (5.15) 

where p{t) is a weak solution of 

div(Vp(t)) = div{p{t)w{t)) in Q, 
Vp{t) • n = on dQ . 

Dividing ()5.15p by p{t), we get 

^i'^) = ~rT'^(*) ^ — rT^P(0 almost everywhere in (0, oo) , (5.16) 
p{t) p{t) 

which after multiplication with Vr for r G H^{i^) and integration gives 

/" 4Tw(t) • Vrdx = - /" -^Vp(t) ■ Vrdx for ah r G F^(17) 
P(^) Jn P{t) 

for almost all t G (0, oo). But this is exactly the weak version (|5.13p of the above auxiliary 
problem {Aux) and we get Vp{t) = Vp{t) for almost every t due to the unique solvability. By a 
redefinition of Vp on a set of measure zero we get the continuity Vp G BCw{[0, T]; L^(f])'^), since 
Vp has this property. Again with a redefinition of v through equation ()5.16p on a set of measure 
zero, we get finally v G BCy,{[0,T]; L'^{nf), where we used the fact p G BUC{[0,T]; L'^{n)) 
with /9 > c > is already known from above. 
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It remains to show that v attains the initial value vq. Due to the above bound of Fo-(/Jv ) in 
Wl{0,T;W^\n)) ^ BUC{[0,T];W^^{n)) we get for arbitrary G C^A^) the convergence 

7 ' 

[ ¥„(^^){0)-iPdx^ [ P,(pov(0)) • V = [ pov(0)-^dx. 
Jn Jn Jq 

By definition of the time-discrete functions we also have F^{pv'^)\t=o = Po-(povo), which together 
with the last convergence yields 

/ povo ■ipdx= I /Oov(O) • ip dx for ah tp e C^^i^) ■ 
Jn Jn 

From Lemma l5. II we get finally v(0) = vq in L^(il)'^. 



5.3 Energy inequality 

Finally we can finish the proof by showing the energy inequality p.l4p . Since v^(t) — )■ v(t) 
in L'^{Q,)'^ and — )• ^p{t) in H^{Q,) for almost every t E (0, oo) (for a subsequence), which 
follows from the strong convergences of and 99^, it holds that 



E^^'it) Etot{<f{t),v{t)) for almost ah t £ (0,oo) . 

Moreover, by lower semicontinuity of norms and almost everywhere convergence of ip^ to (p, the 
inequality 

fOO POO 

lim inf / D'^ {t)T{t) dt > / D{t)T{t)dt 
for all r G W^'^{0,oo) with r > holds, where 

D{t) := / 2r]{ip)\DYf dx + / m((/j)| V/xp . 



Hence, passing to the limit in (jS.Sp . we obtain 

fOO fOO 

£^tot(93o,vo)r(0)+ / Etot{p{t),v{t))T{t)dt> D{t)T{t)dt (5.17) 
Jo Jo 

for all T G Ty^'^(0, 00) with r > 0. With the help of Lemma 12.21 we obtain the energy esti- 
mate IKWi . □ 



A Appendix 

In this last section, we want to discuss existence of weak solutions for a related model discussed 
in Remark 2.2 of Abels, Garcke, Griin [AGGllj . which is derived in detail in [AGGlOj . Although 
this model is not frame indifferent in the standard manner, it can be used to approximate the 
sharp interface model, see [AGGlOj . In addition it has the mathematical interesting term |vp 
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in the equation for the chemical potential and we want to give a short description how the 
methods from this work can be adapted to derive an existence result also in that case. After 
reformulation of the pressure term with the previous notation, this problem is given as 

U.|2 



5t(pv) + div (/9v (g) v) - div {2r]{(f)D\-) + 


= ^Vv3+ 2 Vp 


in 




(A.l) 


divv 


= 


in 


Q, 


(A.2) 


dt(f + V • Vip 


= div (?7i((/?)V/i) 


in 


Q, 


(A.3) 




= ^',{A{^))-^A{^) 


in 


Q, 


(A.4) 




= 


on 


s, 


(A.5) 




= 


on 


s, 


(A.6) 


(v,V')|t=o 


= (vcV^o) 


in 


n. 


(A.7) 



In this case we can show the following existence result. 

Theorem A.l. Let T e (0,oo] and set either I = [0,oo), if T = oo or I = [0,r], if T < oo. 
Let Assumption \3.1\ hold, vq E L'^iS^) o,f^d (po £ H^{i^) with \ipQ\ < 1 almost everywhere and 
-f^ipodx G (—1,1). Then there exists a weak solution (v, of (|A.ip - (|A.7p in the following 
sense: 

V G BC^{I- LliO)) n L2(0, T; HH^Y) , 

^ G BC^{i- hHq)) n LlJL; H^n)) , ^'{^) E L^JL; L\n)) n lIiJI; L\n)) , 

IJieLli^^{L-H\a)) with VfieL\0,T;L\n)) 
and the following equations are satisfied: 

- (P^, dt'^)Q^ + (div(pv v), V')q^ + {2r]i^)Dv, Z)^)q^ 

= (AiVvp, 'iP)q^ + (^Vp, tp^ ^ (A.8) 

for all ip G [C^{^ X (0,T))]'^ wii/i div^ = 0, 

- {if, dtOg^ + (v • Vip, Oq^ = - (m(^)Vp, VC)q^ , (A.9) 
for allC€^C^mT);CH'n)), 

jjL + — —\ + 'kA{lp) = ^o(<^) — Aj4((/9) almost everywhere in Qt and (A. 10) 

(v,V')lt=o = (vcV'o) ■ 

Moreover, 

Etot{p>{t),^{t))+ I 2r]{p)\Dw\^d{x,T)+ I m{p)\V^l\^d{x,T) 

<Etot{^{s)Ms)) (A.12) 
for all t £ [s, oo) and almost all s G [0, oo) holds (including s = 0). 
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The main ideas of the proof of Theorem lA.il are the same as for the previous existence result 
and we will just sketch the necessary adaptations. Again we use an implicit time discretization. 
Leth = jjfoTNe N, G Ll{n), ifk G H\n) with^'{^k) S L'^{n) and pfc = ^{pi+P2)+^{p2- 
pi)(pk be given. We also set ^ = |^ = i(p2 - Pi)- We construct (v, 93,//) = (v^+i, /x^+i) 
as solution of the following non-linear system. 

Find (y,ip,n) with v G H^{nfnLl{n), ip £ V{dE) and p E H^{^), such that 

^ + (div(pfev v), + {2r,{^k)D^r, D^)^ 
= (;uV(/9fc,V)n+ (^Vpfe,v) , (A.13) 

for allt/j G C^A^), 
for ah C e -H'H^^) and 



A(,) - AM V" + "4) + ^^^^^ = - AA(^) a.. i„ a (A.15) 

To show existence of solutions of the time-discrete problem we use again the Leray-Schauder 
principle and define operators Ck,J^k-X^ Y, where 



X = i^H^inf n Ll{Q)j X V{dE) X H^{n) , 

Y = (H^ifif n Llin))' X {H\n)y X L\n) 



For w = (v, if, p) £ X we set 



>Cfc(w) = ( -div(m(v9fc)V/i) + j^pdx 
A{^) + dE{A{^)) 



where 



{Lk{v),ip)= [ 2r]{ipk)Dv : Di/^dx for tp e H^in^ H Ll{n) and 
Jn 

■ / m((/?fe)V/i • V?7 + / pdx ■ / 77 for 77 G //^(il) 
jr2 Jn Jn 



div(m(99fc)V^) + / pdx,!] 
Jn 

Note that due to (/? G T>{dE) it holds G T>{dE) and therefore the last line in Ck{w) lies in 
L^(r2). Furthermore for w = {y,Lp,p) G X we define 



/_pv-^ _ ^.^(^^^ ^ ^ ^^^^ ^ ^^^^ 
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Then w = (v, /u) G X is a weak solution of the time discrete problem (|A.13|) - ()A.15p if and 
only if £fc(w) — J^fc(w) = 0. To get a continuous and even compact operator, we introduce for 
< s < J the following Banach spaces 

X := (^H^iny X Llinj^ X H^-'{n) X H\n) , 

Y :=L^nf X L^{Q) x Wl{n). 

2 

With analogue arguments as above we can show that the inverse : y — t- X is a compact 
operator and that Tk ■ X ^ Y is continuous and maps bounded sets into bounded sets. 

For these operators it is possible to verify the assumptions of the Leray-Schauder principle 
and therefore get a weak solution to the time-discrete problem, which additionally fulfills the 
energy estimate ()4.7p . 

For the remaining part about compactness in time we can show the same statements as 
in Lemma 14.21 and derive then together with the fact P'-^ dx £ L°°(0,oo) from the energy 
estimate the same bounds as in ()5.7p . But in this case the additional term |vp in the line 
(|5.3p for the chemical potential is just bounded in the space Li{J; L^($7)). By applying results 
of Abels and Wilke |AWn7j we get here G LI^^^{J]H'^{Q)) and ^'(99) G LI^^^{J;L'^{Q)). 
The strong convergence ()5.10p can then be derived by using the embedding Ls (0, T; ^^(f^)) ^ 
Li(0, T; H'^{VL)) and the result from interpolation theory V-{Q, T; H'^{n)) D L~(0, T; H\0.)) ^ 
L^{0,T;H2{n)), which shows that is bounded in the latter space and is enough to conclude 
()5.10p with the help of the Lemma of Aubins-Lions. 

To show strong convergence of to v in L^(ilj') in this case, we have to bound in- 
stead of the term (g) in the first model now the term \Y^\'^'V(p^ . From the bound 
of in L2(0,r;L6(Jl)) and L'^{0,T; L'^{n)) we get that jv^p is bounded in L^{0,T; L^{n)) n 
L~(0,T;Li(fi)). With the help of the interpolation results (l5TT]l and (l5T2]l this leads to a 
bound of |v^|2 in lI (0, T; L2(fi)) and together with the bound of Vipf^ £ L°°(0, T; L2(0)) we 

get in this case the boundedness of \v^\^Vip^ in {0,T; (Q)). Anyhow, we have a time 
integrability greater than 1, which is enough to finish the proof. 
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